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ABSTRACT 

We present multiwavelength optical linear polarimetric observations of 69 stars 
toward the young open cluster Be 59. The observations reveal the presence of three 
dust layers located at the distances of ^300, ~500 and ~700 pc. The dust layers 
produce a total polarization Py ~ 5.5 per cent. The mean values of polarization and 
polarization angles due to the dust layers are found to increase systematically with 
distance. We show that polarimetry in combination with the {U — B) — {B — V) colour- 
colour diagram yields a better identification of cluster members. The polarization 
measurements suggest that the polarization due the intra-cluster medium is ~ 2.2 per 
cent. An anomalous reddening law exists for the cluster region, indicating a relatively 
larger grain size than that in the diffuse ISM. The spatial variation of the polarization 
and E(B — V) is found to increase with radial distance from the cluster center, whereas 
the 9v and Xmax are found to decrease with increasing radial distance from the cluster 
center. About 40 per cent of cluster members show the signatures of either intrinsic 
polarization or rotation in their polarization angles. There is an indication that the 
star light of the cluster members might have been depolarized because of non-uniform 
alignment of dust grains in the foreground dust layers and in the intra-cluster medium. 

Key words: Polarization- dust, extinction - open clusters and associations: individ- 
ual: Berkeley 59. 



1 INTRODUCTION 

Interstellar grains are aspherical in nature and, given proper 
conditions, are aligned in space by the magnetic field (Davis 
& Greenstein 1951). The effective extinction cross sections 
of the dust particles are the greatest when the electric vec- 
tor of the incident light is parallel to the long axes of the 
dust particles as projected on the plane of the sky, and the 
least when parallel to the short axes. This differential ex- 
tinction introduces a small degree of linear polarization in 
the transmitted light. 

Studies of polarization due to the interstellar matter 
(ISM) are important as these provide information about the 
properties of the dust associated with the ISM and intra- 
cluster matter as well as help to trace the Galactic mag- 
netic field. As the grains are thought to align due to the 
local magnetic field, the observed polarization vectors map 
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the general geometry of the magnetic field. The observed 
maximum upper limit relation between the degree of the 
polarization and the colour excess E{B — V) is found to be 
Pmax = 9 X E{B-V) (Aannestad & Purcell 1973). The rela- 
tion between Pmax and colour excess, and the variation of P 
with wavelength are interpreted in terms of the grain prop- 
erties and the efficiency of the grain alignment. Therefore, 
polarimetry is a useful technique to investigate the proper- 
ties like maximum polarization P\^^^, the wavelength Xmax 
corresponding to P\^^^ and the orientation of the magnetic 
field in various Galactic locations. 

Polarimetric studies of star-forming regions/young star 
clusters are specially important because physical parame- 
ters such as distance, age, membership and colour excess 
E[B — V) of these regions are known accurately, which 
consequently helps in analyzing the polarimetric data in a 
meaningful way. Strong ultraviolet radiation from 0/B type 
stars in these regions has strong impact on the surrounding 
medium. Dust grains can undergo destruction processes due 
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to direct radiative pressure, grain-grain collisions, sputtering 
or shattering etc. As a result, it is likely that the mean size 
of the dust grains could be smaller than the mean value for 
the diffuse ISM. The stars associated with the star-forming 
regions can help to understand the nature of dust as well as 
the magnetic field of the intra-cluster medium. 

Young star clusters (age < 10 Myr), still embedded in 
the parent molecular clouds, are unique laboratories to un- 
derstand the dust properties as well as the nature of inter- 
action between young star(s) and the surrounding medium. 
Berkeley 59 (02000 = OOh 02m 13s, ($2000= +67° 25' 11"; / 
= 118°. 22, b = 5°. 00) is such a young star cluster associated 
with a heavily obscured gas-dust complex of the Cepheus 
OB4 association. The cluster Be 59, located at the center 
of the Sharpless region S171 contains nine 07— B3 stars (cf. 
Pandey et al. 2008; hereafter P08) at a distance of 1.00±0.05 
kpc and has E{B — y) ~ 1.4 — 1.8 mag. The extent of the 
cluster was found to be 2.9 pc (P08). 

As a part of an ongoing project to understand the dust 
characteristics in star-forming regions and to map the struc- 
ture of magnetic field at diverse environments of the Milky 
Way Galaxy, we have carried out broad-band optical polari- 
metric observations around the cluster Be 59. In Section 2 
we present the observations and data reduction. Results are 
presented in Section 3. The dust properties and the spatial 
variation of E[B — V), Py, 9v and \max are in Sections 4 and 
5 respectively. Finally, we conclude our results in Section 6. 



tion using the IRAPy package. The ratio R{a) is computed 
using the following relation 

Ri^)= 'j:t] ^ =Pcos{2e-4a), (1) 

where P is the fraction of the total light in the linearly po- 
larized condition and 9 is the position angle of the plane of 
polarization. Here a is the position of the fast axis of the half 
wave plate (HWP) at 0°, 22.5°, 45° and 67.5° correspond- 
ing to the four normalized Stokes parameters respectively, 
q [R(0°)], u [R(22.5°)], gi [R(45°)] and m [R(67.5°)]. The 
detailed procedures used to estimate the polarization and 
position angles for the program stars are given by Eswara- 
iah et al. (2011) (hereafter Ell) and references therein. 

The instrumental polarization of AIMPOL on the 104- 
cm Sampurnanand Telescope has been monitored since 2004 
on various observing nights and found to be less than 0.1 per 
cent in different bands (Ell and references therein). All the 
measurements were corrected for both the null polarization 
(~ 0.1 per cent) which is independent of the pass-bands and 
the zero-point polarization angle by observing several unpo- 
larized and polarized standard stars from Schmidt, Elston & 
Lupie (1992) (hereafter S92). The results for the standard 
stars are given in Table [T] The first column lists the star 
name with the date of observation or reference in paren- 
thesis. The next consecutive columns are the polarization 
in percent [P (per cent)] and polarization angle in degree 
[0{°)] measured in BV{RI)c pass-bands. Entries with S92 
in parenthesis are taken from S92. The present results for 
the polarized standard stars are in good agreement, within 
the observational errors, with those given by S92. 



2 OBSERVATIONS AND DATA REDUCTION 

Polarimetric observations were carried out on seven nights 
(2009 November 23, 24, 25 and 2009 December 24, 26, 
27 and 28), using the ARIES Imaging Polarimeter (AIM- 
POL; Rautela, Joshi & Pandey, 2004) mounted at the 
Cassegrain focus of the 104-cm Sampurnanand telescope of 
the Aryabhatta Research Institute of observational sciencES 
(ARIES), Nainital, India. The observations were carried out 
in the B, V, Rc and 7c (As^j,j:=0.440/im, Av;.J:J^=0.530/im, 
Aflc^jj:=0.670/im and A/^^j, =0.800/im) photometric bands 
using fraction (370x370 pixel^) of the TK 1024 x 1024 
pixel^ CCD camera. The AIMPOL consists of a half-wave 
plate modulator and a WoUaston prism beam-splitter. The 
WoUaston prism analyzer is placed at the backend of the 
telescope beam path in order to produce ordinary and ex- 
traordinary beams in slightly different directions separated 
by 28 pixels along the North-South direction on the sky 
plane. A focal reducer (85 mm, f/1.8) is placed between the 
WoUaston prism and the CCD camera. Each pixel of the 
CCD corresponds to 1.73 arcsec and the field of view is '--^ 8 
arcmin diameter on the sky. The FWHM of the stellar im- 
age varied from 2 to 3 pixels. The read-out noise and gain of 
the CCD are 7.0 e"^ and 11.98 e"VADU respectively. Since 
AIMPOL does not have a grid, we manually checked for any 
overlap of ordinary and extraordinary images of the sources. 
Fluxes of ordinary (7o) and extra-ordinary (7e) beams for all 
the observed sources with good signal-to-noise ratio were ex- 
tracted by standard aperture photometry after bias subtrac- 



3 RESULTS 

Table [5] lists the polarization measurements for 37 stars us- 
ing BV{RI)c bands, whereas Table|3]lists the results for 32 
stars using the V{RI)c bands. The star identification num- 
bers (column 1) are taken from P08. The right ascension, 
declination and photometric visual magnitudes, also from 
P08, are listed in the second, third and fourth columns of 
Tables [2] and O respectively. The next consecutive columns 
correspond to the polarization, polarization angle and their 
associated errors in B, V, {R, I)c bands and V, {R, I)c bands 
respectively. The given polarization angles are in the equato- 
rial coordinate system measured from the North increasing 
towards the East. Tables [2] and |3] reveal that the maximum 
linear polarization towards the cluster region is ~ 8 per cent. 
Such a high amount of polarization is not often found to- 
wards the star clusters, with only a few exceptions, e.g., Tr 
27 (Feinstein et al. 2000), and M17 (Schultz et al. 1981). 

In Fig.[T]all the observed 69 stars are marked with white 
circles on the DSS II 7?-band image. The sky projection of 
l/-band polarization vectors are overlaid. The length of each 
polarization vector is proportional to the degree of polariza- 
tion. The dash-dotted line represents the orientation of the 
projection of the Galactic plane (GP) at fe=5.03°, which cor- 
responds to the position angle of 86° . 

IRAF is distributed by National Optical Astronomical Obser- 
vatories, USA. 



Table 1. Observed polarized standard stars from S92 



star name Pb ± e 6b ± e Py ± e By ±e Puc ± e ^jjc i £ -P/c i e Ojc =t e 
(date of observations) 

(reference) (per cent) (°) (per cent) (°) (per cent) (°) (per cent) (°) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



HD 19820 (23, Nov 09) 
HD 19820 (24, Nov 09) 
HD 19820 (25, Nov 09) 
HD 19820 (23, Dec 09) 
HD 19820 (24, Dec 09) 
HD 19820 (27, Dec 09) 
HD 19820 (28, Dec 09) 
HD 19820 (S92) 



4.49 ± 0.11 
4.61 ± 0.11 
4.72 ± 0.11 
4.70 d= 0.10 
4.68 ± 0.11 
4.57 ± 0.11 
4.72 ± 0.09 
4.70 ± 0.04 



114.9 ± 0.7 
114.7 ± 0.7 

115.6 ± 0.7 
115.0 ± 0.6 
116.4 ± 0.7 
116.4 ± 0.6 
115.4 ± 0.5 

115.7 ± 0.2 



4.89 ± 0.09 

4.79 ± 0.08 
4.94 ± 0.08 

4.90 ± 0.09 
4.69 ± 0.09 

4.80 ± 0.08 
4.79 ± 0.03 



114.2 ± 0.5 
115.1 ± 0.5 
114.9 ± 0.4 

114.6 ± 0.5 
115.0 ± 0.5 

114.8 ± 0.5 

114.9 ± 0.2 



4.49 ± 0.09 
4.51 ± 0.07 
4.61 ± 0.07 

4.61 ± 0.07 

4.62 ± 0.07 

4.63 ± 0.07 
4.51 ± 0.07 
4.53 ± 0.03 



115.5 ± 0.6 

114.6 ± 0.4 

114.4 ± 0.4 
114.3 ± 0.4 
114.8 ± 0.4 
114.6 ± 0.4 
114.2 ± 0.4 

114.5 ± 0.2 



4.06 ± 0.16 
3.97 ± 0.09 
3.86 ± 0.10 
4.06 ± 0.08 
4.06 ± 0.08 
4.01 ± 0.08 
3.92 ± 0.08 
4.08 ±0.02 



115.3 ± 1.0 
115.7 ± 0.6 

113.4 ± 0.7 
114.0 ± 0.5 
114.3 ± 0.5 
114.3 ± 0.5 
115.2 ± 0.5 

114.5 ± 0.2 



HD 25443 (23, Nov 09) 
HD 25443 (24, Nov 09) 
HD 25443 (25, Nov 09) 
HD 25443 (23, Dec 09) 
HD 25443 (24, Dec 09) 
HD 25443 (27, Dec 09) 
HD 25443 (S92) 



5.19 ± 0.09 
5.21 ± 0.09 

5.11 ± 0.09 

5.12 ± 0.08 
5.12 ± 0.09 
5.05 ± 0.08 
5.23 ± 0.09 



134.6 ± 0.5 
134.1 ± 0.5 
134.5 ± 0.5 

135.7 ± 0.5 
134.0 ± 0.5 

134.8 ± 0.5 
134.3 ± 0.5 



5.04 ± 0.07 

5.11 ± 0.07 
5.25 ± 0.09 

5.12 ± 0.08 
5.27 ± 0.09 
5.17 ± 0.07 

5.13 ± 0.06 



136.0 ± 0.4 

136.2 ± 0.4 

134.3 ± 0.5 
135.5 ± 0.4 

134.8 ± 0.5 

135.9 ± 0.4 
134.2 ± 0.3 



5.01 ± 0.06 
5.13 ± 0.06 
5.10 ± 0.08 
4.88 ± 0.08 
5.00 ± 0.08 
5.04 ± 0.07 
4.73 ± 0.05 



134.5 ± 0.4 
133.8 ± 0.4 

133.4 ± 0.4 

134.5 ± 0.5 

134.5 ± 0.4 
134.0 ± 0.4 

133.6 ± 0.3 



4.19 ± 0.09 
4.29 ± 0.09 
4.37 ± 0.10 
4.23 ± 0.08 
4.19 ± 0.07 
4.13 ± 0.07 
4.25 ± 0.04 



134.8 ± 0.6 

134.4 ± 0.6 
132.6 ± 0.7 

134.9 ± 0.5 
136.0 ± 0.5 

134.5 ± 0.5 
134.2 ± 0.3 



BD+64 106 (23, Nov 09) 
BD+64 106(23, Dec 09) 
BD+64 106 (S92) 


5.49 ± 0.17 
5.37 ± 0.15 
5.51 ± 0.09 


98.0 ± 0.9 
97.9 ± 0.8 
97.2 ± 0.5 


6.09 ± 0.13 
5.57 ± 0.11 
5.69 ± 0.04 


99.0 ± 0.6 
96.8 ± 0.6 
96.6 ± 0.2 


5.41 ± 0.11 
5.63 ± 0.09 
5.15 ± 0.10 


96.1 ± 0.6 

97.6 ± 0.5 

96.7 ± 0.5 


4.50 ± 0.14 
4.63 ± 0.10 
4.70 ± 0.05 


96.6 ± 0.9 
97.5 ± 0.6 
96.9 ± 0.3 


HD 204827 (25, Nov 09) 
HD 204827 (S92) 


5.67 ± 0.10 
5.65 ± 0.02 


58.3 ± 0.5 
58.2 ± 0.1 


5.53 ± 0.08 
5.32 ± 0.01 


58.8 ± 0.4 
58.7 ± 0.1 


5.02 ± 0.08 
4.89 ± 0.03 


59.2 ± 0.4 
59.1 ± 0.2 


4.14 ± 0.10 
4.19 ± 0.03 


59.9 ± 0.7 
59.9 ± 0.2 


BD+59 389 (23, Dec 09) 
BD+59 389 (28, Dec 09) 
BD+59 389 (S92) 


6.35 ± 0.13 
6.43 ± 0.13 
6.34 ± 0.04 


97.8 ± 0.6 

97.9 ± 0.6 
98.1 ± 0.2 


6.73 ± 0.09 
6.82 ± 0.09 
6.70 ± 0.01 


97.8 ± 0.4 
97.5 ± 0.4 
98.1 ± 0.1 


6.48 ± 0.08 
6.47 ± 0.08 
6.43 ± 0.02 


97.6 ± 0.3 

97.7 ± 0.4 
98.1 ± 0.1 


5.66 ± 0.06 
5.61 ± 0.07 
5.80 ± 0.02 


97.9 ± 0.3 
97.9 ± 0.4 
98.3 ± 0.1 


HD 236633 (28, Dec 09) 
HD 236633 (S92) 


6.06 ± 0.11 
5.53 ± 0.04 


90.4 ± 0.5 

92.5 ± 0.2 


5.65 ± 0.09 
5.49 ± 0.02 


91.3 ± 0.4 
93.8 ± 0.1 


5.34 ± 0.09 
5.38 ± 0.03 


91.0 ± 0.5 
93.0 ± 0.2 


4.69 ± 0.09 
4.80 ± 0.04 


90.5 ± 0.6 
93.1 ± 0.2 
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Figure 1. The stellar polarization vectors are superimposed on a 21' X 19' R-band DSS II image of the field containing Be 59. The 
length of the polarization vector is proportional to Py ■ A vector with a polarization of 6 per cent is drawn for reference. The dash-dotted 
line is the Galactic plane at b=5.03°. The plus symbols represents the cluster center RA=00''02'"10*.4 and Dec=67°25'10". The cluster 
radius is shown with a big circle of ~10 arcmin (P08). The stars observed are identified using the star identification numbers from P08. 
North is at the top and east is to the left. 



Table 2. Observed BV(ii/)c polarization and polarization angles for 37 stars towards Be 59 



Star ID* R.A (°) DEC (°) V (mag)* Ps ± e ± e Pv±e dy ± e Prc±<^ f Ac ± e P/c ± e 9i^±e 

(2000J) (2000J) (per cent) (°) (per cent) (°) (per cent) (°) (per cent) (°) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 
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* Star IDs and V magnitudes, barring 16A, have been taJjen from P08 



Table 3. Observed V, {RI)c polarization and polarization angles for 32 stars towards Be 59 
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3.1 Distribution of Pv and 6v- 

A careful inspection of Fig. [1] reveals two groups of stars 
characterized by different degree and direction of polariza- 
tion. The first group with relatively small degree of polariza- 
tion (~ 2 per cent) and with the orientation nearly parallel 
to the Galactic disk 86°) may be composed of foreground 
stars. The stars of this group are randomly distributed on 
the plane of the sky. The second group, whose degree of po- 
larization is significantly higher than the first group and the 
alignment of polarization vectors (~ 102°) is significantly 
deviated from the GP, may be composed of cluster mem- 
bers. Majority of the stars of the second group are located 
within the cluster region (~ 10 arcmin). 

Figure [2] displays the distribution of polarization Py vs 
polarization angle 9v, which clearly segregates the above 
mentioned two groups. The mean Py and Oy for the first 
group is found to be 2.47 ±0.46 per cent and 81°±9°, re- 
spectively. The mean polarization angle is nearly aligned 
with the GP, indicating the homogeneity in the local mag- 
netic field of the Galaxy towards the direction of Be 59. 
The stars belonging to the second group are more polarized 
(Pv > 4.0 per cent) and having 9v ^ 90°. The mean po- 
larization angle of the second group (103° ±5°) significantly 
deviates from the GP. The magnetic field associated with 
the parental molecular cloud may have been perturbed dur- 
ing the cloud collapse or due to the strong stellar winds or 
supernova explosions (e.g., Waldhausen, Martinez & Fein- 
stein 1999). The stars belonging to the second group could 
be either the cluster members or background stars. The sec- 
ond group of stars show large dispersion in polarization but 
a very small dispersion in polarization angle. The large dis- 
persion in Pv (4 — 8.7 per cent) could be attributed to the 
differential reddening within the cluster. Similar type of seg- 
regated groups have already been reported in a few cases, 
namely Trumpler 27 (Feinstein et al. 2000), Hogg 22 & NGC 
6204 (Martinez et al. 2004), NGC 5749 (Vergne et al. 2007) 
and NGC 6250 (Feinstein et al. 2008). 

The distribution of Pv for all the observed 69 stars is 
shown in the upper panel of Fig. [S] which clearly reveals two 
separate distributions for field and probable cluster mem- 
bers. The distribution for field stars peaks at ~ 2 per cent, 
whereas that for the probable members peaks at 5.5 per cent 
with an extended tail towards higher polarization, which 
could be either due to highly extincted probable members, 
background stars or due to the presence of different popula- 
tions of dust grains with different polarizing properties. The 
lower panel of Fig. [3] shows the distribution of 9v, which re- 
veals that the distribution of probable cluster members lies 
in the range of 95° — 115° with a peak at ~105°. 

To have a better understanding of the nature of the 
dust component and the magnetic field associated with the 
foreground and intra-cluster medium, it is essential to find 
out the members associated with the cluster versus the stars 
located in the foreground/background of the cluster. In our 
previous study (Ell) we have shown that polarimetry in 
combination with the {U — B) — {B — V) two-colour dia- 
gram (TCD) can yield a better identification of probable 
members than photometry alone. In the ensuing section we 
will discuss the determination of membership using the po- 
larization properties in combination with (U — B) — {B — V) 
TCD. 



120 F 




P. (%) 

Figure 2. Polarization angle versus polarization in V-band for 
69 stars towards Be 59. Dashed lines are drawn to show the two 
clearly separated groupings among the observed sample. 
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Figure 3. Histograms for Pv (upper panel) and dv (lower panel) 
for 69 stars towards Be 59. Thick line drawn in the lower panel 
represents the projection of the Galactic plane at 86°. 
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open and filled circles. The late type stars are shown with filled diamond symbols. 



3.2 Member identification 

[U ~ B) ~{B ~ V) TCD is one of the useful tools to iden- 
tify probable members of a cluster. It is expected that all 
member stars have E[B — V) values comparable to the mean 
E[B — V) value of the cluster since the cluster stars have 
formed out of the same molecular cloud and consequently 
have the same distance and age. In comparison, the field 
population is expected to be less or highly extincted de- 
pending on whether they are foreground or background to 
the cluster. Fig. g] shows the {U - B) - {B - V) TCD for 
only 63 stars as the U — B colours are not available for 6 
stars. In Fig. [J] the zero-age-main-sequence (ZAMS) from 
Schmidt-Kaler (1982) is shifted along a normal reddening 
vector having a slope of E{U - B)/E{B ~V) ^ 0.72. The 



TCD shows a variable reddening in the cluster region with 
E{B - V)mi„~1.4 mag and E{B - V)max~l.^ mag. The 
cluster members (shown by open circles) seem to have spec- 
tral types earlier than AQ. The TCD manifests the presence 
of a foreground population shown by filled circles. The fore- 
ground population reddened by E[B — V)=QA mag is found 
to be located at ~ 470 pc (P08). The other field star (late 
type) population is shown by filled diamonds. 

The polarimetric observations can be used as a tool to 
identify member stars in a Galactic open cluster, particu- 
larly, when the field stars have colours similar to those of 
cluster members (Martinez et al. 2004; Vergne et al. 2007; 
Feinstein et al. 2008; Vergne et al. 2010; Orsatti et al. 2010; 
Ell). The individual Stokes parameters of the polarization 
vector of the l/-band, Pv, given by Qv = Pv cos{26v) and 
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Uv = Pv sm{26v) are estimated for all the observed stars 
toward Be 59 and presented on a Uv versus Qv plot, known 
as the Stokes plane, in Fig. (5] 

The measured degree of polarization of a star depends 
on the cumulative amount of aligned dust grains that lie 
along the line of sight, and hence the degree of polarization 
would be similar, lower or higher depending on whether it is 
member, foreground or background to the cluster. Likewise, 
the position angles of the cluster members would be simi- 
lar, but different for foreground or background field stars as 
light from them could have contributions from different or 
additional dust components. Hence the cluster members are 
expected to show a grouping in the Uv — Qv plane, while 
non-members are expected to show a scattered distribution. 
Therefore Uv — Qv plot could be a useful tool to identify 



the members of a cluster. The stars with intrinsic polariza- 
tion (due to asymmetric distribution of matter around young 
stellar objects (YSOs)) and/or rotation in their polarization 
angles may also create scattered distribution in the Uv — Qv 
plane. Nebulous background in case of star-forming regions 
would also create small intrinsic polarization and hence pos- 
sibly the scattered distribution. 

Fig. [5] shows two prominent groupings around Uv ~ 
—2.5 and Qv ~ —5.5 (first group) and Uv ^ 1.0 and Qv ~ 
—2.0 (second group). The grouping at Uv ~ —2.5 and Qv ^ 
—5.5 should be due to the cluster members. Fifty percent (13 
out of 27) of the field stars, identified on the basis of {U — 
B) - {B - V) TCD (Fig. HI, can be noticed mainly around 
Uv ~ 1.0 and Qv ~ —2.0. The remainder of the probable 
field stars show a scattered distribution in Uv ~ Qv plot; 
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however some of the probable field stars are found to mingle 
with the probable cluster members. To further elucidate the 
membership we plot a box with dashed line in Uv — Qv plot 
having boundaries of mean Py ± c^p^, (6.22±1.21 per cent) 
and mean By ± cre^ (102± 6°) obtained using the probable 
member stars (open circles) identified in Fig [J] The stars 
shown with open circles and lying within the Icr box of the 
mean Pv and 6v could be probable members of the cluster. 
It is apparent from Fig. H i.e., the {U - B) - {B - V) TCD 
that the majority of the stars located within the Icr box 
follow the general reddening of the cluster region, hence are 
probable members of the cluster. 

Stars # 59, 120, 689, 828 and 1213 are located signif- 
icantly away from the Icr box. The stars #59 (P\/=4.5 per 
cent, 0^=90°) and #120 (Pv=3.0 per cent, 6lv = 107°), even 
though located spatially within the cluster region (see Fig.[T]) 
and with photometric colours consistent with membership, 
their Py value or 6v value is not comparable to the cluster 
region (see Fig. [5]), so are considered as non-members. The 
stars # 689, 828 and 1213 are located outside the cluster 
region. These stars show relatively large polarization (~ 6.6 
to 8.7 per cent) but their By values (~ 88° to 91°) are sig- 
nificantly different from those of probable cluster members 
(Tables[2]and[21) and are comparable to the GP. Hence these 
are considered as non-members. The stars # 23, 195, 247, 
324, 391, 418, 441, 433 and 461 are distributed outside but 
near the boundary of la box. Of these 9 stars (barring # 23, 
195 and 324) have relatively higher value of e (cf. Fig. IIOI 
Sec[4| indicating rotation in their polarization angles. These 
stars have E{B — V) in the range of 1.4—1.8 mag (cf. Fig. 
[l]), so could be members of the cluster. Three stars, namely 
# 20, 62 and 367 are also located near the boundary of the 
Icr box. The Pv values of these stars are in the range of 4.3 
per cent to 4.6 per cent which is the lower limit of Pv value 
for the cluster member stars (see Tables [2] and [3| . The 6v 
values range from 105° — 108°. These could also be members 
of the cluster. 

There seems to be a less prominent grouping around 

Uv 1.0 and Qv ~ -5.0 consisting of stars # 202, 742, 

936, 1067, 1273, 1305 and 1684 (cf. Fig.[5|. We refer to this 
group as the third group. Interestingly, all these stars (ex- 
cept 202) are located spatially at the same region towards 
the Northern part of the cluster Be 59. This group of stars 
is located near the edge of Icr box. The mean value of Pv 
(5.68±0.53 per cent) is comparable to the Pv values of the 
cluster region. However, the mean 6v (=95° ±2°) is signifi- 
cantly different from that in the cluster region. These stars 
are considered as field stars. 

The colours of stars # 1028, 1135 and 1451 are com- 
parable to those of the cluster members and lie in the 
E[B — V)~1.4 to 1.8 mag range (presuming that these stars 
have spectral types earlier than ^0). Star # 1135 has Pv 
(8.12±0.45 per cent) higher than those for cluster members, 
however 6v (96.6±1.5°) is smaller than the value for clus- 
ter stars. This star is located outside the Icr box and also 
located outside the estimated boundary of the cluster (Fig. 
[U, hence this star is considered as a non- member. Stars # 
1028 and 1451 have Pv and 9v values both comparable to 
those of the cluster members but located outer side bound- 
ary of the cluster, hence the membership of these stars is 
uncertain. 

Star # 196 is located near the solar neighborhood in the 



Qv and Uv diagram. Its colours are consistent with those 
of early type members of the cluster, which yields E{B — 
V) = 1A mag, but its location in the Qv — Uv diagram (Fig. 
[5]) manifests that it should be a field star, hence this star is 
considered as a field star for further discussion. The star # 
1106 is located near the Icr box boundary (cf. Fig.[S)). On the 
basis of its {U — B) ~ {B ~V) colours and its location outside 
the cluster boundary, we consider it as a non-member. 

The stars # 382, 503 and 988 are well separated from 
the cluster probable members in the {U — B)-{B — V) TCD 
(cf Fig. |4]), however they lie within la box in the Qv — Uv 
plot (Fig.[S}. Since stars # 382 and 503 lie within the bound- 
ary of the cluster and Pv and 9v values are comparable to 
the cluster stars, these are considered as probable members. 
Star # 988 lies outside the boundary of the cluster, its mem- 
bership is uncertain. 

The probable members of the cluster identified using 
Uv — Qv and colour-colour diagrams are given in Tabled 
The member, probable member and non-member stars are 
represented with "M", "PM" and "NM" respectively. Stars 
with uncertainty in their membership determination are in- 
dicated with a "?" symbol. 

Since cluster members seem to have spectral types ear- 
lier than AO (see Fig. |4]), the reddening E{B - V) for 
the member stars has been estimated using the Q-method 
(Johnson & Morgan 1953). As seen in Fig. H) all the field 
stars have spectral types later than AO. The reddening 
E{B — V) for field stars is estimated visually using the slide- 
fit method of ZAMS along the reddening vector. The esti- 
mated values of E{B — V) are given in Table |4] The colours 
of stars # 828 and 1213 suggest spectral types earlier than 
AO. The colours of star # 689 suggest two reddening values, 
{E{B — V)=0.70 mag for a spectral type later than ^40, or 
E{B — V) = 1.44 mag for a spectral type earlier than AO). 
Star # 1135 has Pv'=8.72 per cent, which is comparable to 
the cluster members, however its 9v value (96.6°) is different 
from the cluster region. Its colours also suggest two values 
of E{B — V) (1.17 mag for a spectral type later than AO, 
or 1.42 mag for a spectral type earlier than AO). Similarly, 
colours for stars # 1028, 1106, and 1451 have two possible 
values of reddening. These are mentioned in Tabled] 

The colours of the third group of stars (# 742, 936, 
1067, 1273, 1305 and 1684) suggest that their E{B - V) 
values are in the range of ~ 0.95—1.30 mag. The location 
of the star # 936 in the {U - B) - {B - V) TCD suggests 
an O/B spectral type. The mean values of Pv and 6v for 
these stars are 5.69 ± 0.58 per cent and 96° ±2° respectively. 
This group of stars may lie between the foreground stars (~ 
470 pc) and the cluster. Assuming a mean E{B — V)=1.20 
mag, the Vo/{V — I)o CMD indicates that these stars are 
distributed at a distance of ~ 700 pc. 

Fig. 1^ shows the distribution of Pv (left panels) and 0v 
(right panels) for the identified foreground stars (first group 
and third group of stars) and cluster members. The mean 
and standard deviation of Pv and 9v of these three groups 
of stars are given in a Tableland also shown in Fig.|6] The 
mean values of Qv and Uv are also listed in Table [5] Fig. 
[5] shows the mean degree of polarization, polarization angle 
as well as the deviation of the mean polarization angle from 
the GP increases systematically with increasing distance. 
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Table 4. The E{B — V) values estimated using the {U — B) — (B — V) TCD. The membership information is also mentioned against 
each star ID 



stars with BV{RI)c pass 


-band data 


stars with V{Rr)c pass- 


-band data 


Star ID * 


E(B — V") 


membership 


Star ID * 


E(B ~V) 


membership 


C-) 




(3j 






V°) 


3 


1.41 


M 


47 


0.43+ , 1.20 


NM 


15 


1.38 


M 


59 


0.63+, 1.13+ 


NM 


16 




NM 


62 


1.35 


M 


16A 




M 


64 


1.55 


M 


20 


0.48t, 1.40t 


M 


92 


1.64 


M 


23 


1.80 


M 


105 


1.53 


M 


66 


1.79 


M 


120 


0.70+, 1.06+ 


NM 


69 


1.33 


M 


133 


1.63 


M 


114 


1.52 


M 


195 


1.47 


M 


130 


1.87 


M 


196 


0.57+ 


NM 


157 


1.54 


M 


198 


0.40+ 


NM 


169 


1.55 


M 


202 




NM 


186 




M 


215 


0.63+ 


NM 


234 


1.41 


M 


240 




NM 


239 


0.33t, 0.65t 


NM 


317 


1.81+ 


NM 


247 


1.67 


M 


433 


1.74 


M 


307 


0.95t 


NM 


503 


1.30 


PM 


310 


1.58 


M 


594 


0.35+, 0.51 


NM 


316 


0.65t 


NM 


689 


0.70+ , 1.44 


NM 


324 


1.61 


M 


742 


0.95t, 1.61 


NM 


367 




M 


828 


1.54 


NM 


368 


0.55+ 


NM 


1028 


1.14+, 1.45 


7 


382 


1.26 


PM 


1106 


1.17 


NM 


391 


1.65 


M 


1135 


1.17+ , 1.42 


NM 


418 


1.41 


M 


1190 


0.46+ 


NM 


441 


1.53 


M 


1273 


0.56+, 1.16+ 


NM 


454 


1.42 


M 


1305 


0.95+, 1.17+ 


NM 


461 


1.36 


M 


1312 


1.34+ 


NM 


486 




NM 


1451 


1.08+, 1.40 


7 


936 


1.50 


7 


1513 


0.30+ 


NM 


988 


1.05t 


7 


1523 


0.63+ 


NM 


1067 


0.55t, 1.27 


NM 


1569 


0.45+ 


NM 


1154 


0.85t 


NM 








1213 


1.63 


NM 








1446 


0.45+ 


NM 








1580 




NM 








1684 


0.95+ NM 











* : From P08 

+ : E{B — V) values were obtained using slide and fit method 
M; Members 
NM: Non-mcmbcrs 
PM: Probable members 
?: Stars with uncertainty in their membership 

Table 5. The estimated mean values of Py, 6y, Qv and Uy for the foreground stars, the stars in 3rd group and the cluster members. 
The distance information is also mentioned 



















Py (per cent) 


ov n 


Qv 


Uv 


No. of stars 


Distance (pc) 


First group 


2.44 ± 0.45 


82±10 


Q/=-2.35 


Uf=om 


14 


470 


Third group 


5.55 ± 0.41 


95±2 


Q3=-5.43 


Us = -1.04 


6 


700 


Member stars 


6.34 ± 1.05 


105±4 


Q,„=-5.45 


Um -3.24 


29 


1000 
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3.3 Dust distribution 
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Figure 6. Histograms of Py (left panels) and By (right panels) 
for the foreground stars, stars in 3rd group and the cluster mem- 
bers. For comparison, the GP is also drawn with thick line at 86°. 
The mean and standard deviation of Py and 9y values for each 
group are also mentioned. 
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Figure 7. Uy versus Qy diagram using the stars with known 
membership. The stars near the Sun, 3rd group of stars and the 
cluster members are shown with filled circles, filled squares and 
open circles respectively. The mean Stokes parameters for these 
three groups of stars are also shown with bigger open circle, open 
square and filled circle symbols. The filled star symbol denotes 
the position of Sun at Qy=0 and Uy=0 



The stokes plane can be used effectively not only to deter- 
mine the membership but also to delineate the variations in 
the interstellar environments, e.g., the distribution of dust 
layers, the role of dust layers in polarization, the associated 
magnetic field orientation etc. The vector that connects two 
points in the Stokes plane represents the amount of polariza- 
tion, while any change in the direction of vectors is related to 
change in the polarization angle as seen by us toward a par- 
ticular line of sight. If the dust grains are oriented/aligned 
uniformly (i.e., uniform magnetic field orientation) then the 
degree of polarization is expected to increase with distance 
but the direction of polarization (polarization angle) should 
remain the same and hence the Stokes vector should not 
change its direction with increasing distance. For example, 
in the case of NGC1893 (Ell) the degree of polarization 
was found to increase with distance whereas the direction of 
polarization remains almost constant (cf. their Fig. 5). 

To understand the dust distribution towards Be 59, we 
compare the polarization measurements of foreground stars, 
the stars of third group, and the cluster members as shown 
in Fig. [T] The mean Stokes parameters of these three groups 
of stars are also marked with large open circle, open square 
and filled circle respectively and are connected with dashed 
lines. 

The individual polarization properties of the dust layers 
have been estimated by subtracting their foreground contri- 
bution and the same are given in Table |6l The polariza- 
tion measurements of first group of foreground stars (mean 
Pv'=2.44±0.45 per cent and mean 6lv=82±10°) located at 
^ 470 pc indicate the presence of a dust layer at < 470 pc. 
The orientation of magnetic field of this layer is found to 
be comparable to the GP. The third group of stars, lying 
between the first group and the cluster Be 59, show differ- 
ent polarization measurements (Pv=5.55±0.41 per cent and 
6v=95°±2°) from those of first group or the cluster mem- 
bers. This fact manifests that there must be another dust 
layer at the distance >500 pc that polarizes the star light of 
the second group of stars by ~ 3.5 per cent (see Table [6]). 
The dust grains in this dust layer are found to be aligned sig- 
nificantly differently ('^ 104° , see Table [6]) from that in the 
first dust layer (~ 82°). The polarization measurements for 
cluster members (Pv=6.34±1.05 per cent and 6y= 105±4°) 
also indicate different polarization properties of intra-cluster 
dust. Table [6] suggests that intra-cluster medium also polar- 
izes the cluster members by ~ 2.20 per cent. The dust grains 
of the intra-cluster material are aligned significantly differ- 
ently (~135°) from those in the two foreground dust layers. 
The polarization angles found to increase with increase in 
distance from the Sun (cf. Table |6] and Fig. (Tj . This sys- 
tematic change in the alignment of dust grains may cause 
depolarization (less polarization efficiency). This issue will 
be discussed in the ensuing section. 

To further study the dust distribution for distance < 
500 pc towards the direction of Be 59, we select stars with 
polarization measurements from the catalogue by Heiles 
(2000) and having Hipparcos parallax measurements (van 
Leeuwen, 2007) in a 10° radius around Be 59. We selected 
stars by applying the following criteria: (a) Py> 0.1 per cent, 
(b) ratio of parallax error to the parallax, i.e., (Jt^^/tih ^ 
0.5 and (c) stars without having any emission features or 
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Table 6. The estimated net mean values of Py, By , Qy and Uy due to the dust layers in front of first and third group of stars as well 
as due to the intra-cluster medium. 













Dust layers 


Py (per cent) 


dv n 


Qv 


Uy 


In front of first group of stars 


2.44 


82 


Qf=-2.35 


Uf=0.65 


In front of third group of stars 


3.51 


104 


Qs -Q/=-3.08 


U3 - Uf= -1.69 


Intra-cluster medium 


2.20 


135 


Qm - Q3=-0.02 


Um - C/3=-2.20 
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Figure 8. Distance versus polarization (upper panel) for the stars 
in the direction to Be 59. Polarization and distance information 
for the stars located at distance < 500 pc is obtained from Heiles 
(2000) and van Leeuwen (2007) respectively. Polarization versus 
polarization angle (lower panel). 



photometric variability (with the help of SIMBAD). Fig- 
ure [S] shows polarization versus distance (upper panel) and 
polarization versus polarization angle (lower panel) for the 
stars studied in the present work (same symbols as in Fig. 
O as well as the stars (shown with dots) from Heiles (2000). 
Figure [8] indicates a sudden increase in the polarization at 
~300 pc, ~500 pc and ~700 pc, which suggests the pres- 
ence of three dust layers at ~ 300 pc, ~500 pc and ~ 700 
pc towards Be 59. 

As shown in the lower panel, the polarization angles 





Pm.»(VRl) (%) 



0.4 0,6 
K... (VRI) (^m) 



Figure 9. Left panel: Pmax computed using four passbands (us- 
ing K.=1.6&\max) versus the same parameter but using three 
passband data (using K=1.15). Right panel: Same as left panel 
but for \max- Only 37 stars are used which have four {BV {RI)c) 
pass bands data (cf. Table [Tjl. 



of the Heiles stars (dots) are distributed randomly, which 
indicates that the magnetic field orientation in the nearby 
but in an extended region towards the direction of Be 59 is 
not as organized as compared with that in the cluster region. 
Hence the magnetic field in the intra-cluster medium seems 
to be more confined. 

Neckel & Klare (1980) have studied the reddening dis- 
tribution in the GP with |6|< 7.6° using the extinction and 
distances computed for individual stars. The Ay map to- 
wards the direction of Be 59 by Neckel & Klare (1980) [see 
their Figure 6a, 4 (115/3)] shows an increase in Ay by ~ 
0.6 mag at the distance of 300 pc, indicating the pres- 
ence of a dust layer at this distance. For a normal redden- 
ing law. Ay ~ 0.6 corresponds to E{B — V) ~ 0.20 mag. 
This value of E{B — V) yields a polarization of ~ 1 per cent 
(P=5xiJ(i? — y)) which is in accordance with the dust layer 
at ~ 300 pc as shown in Fig. [8] At ~ 800 pc the Ay further 
increases and reaches ~ 1 mag. The Ay has a steep rise after 
800 pc and at a distance of 1 kpc the Ay is ~ 3 mag corre- 
sponding to the E{B — V') of ~ 1 mag, which is consistent 
with the cluster's foreground reddening. The present polari- 
metric results are therefore consistent with the reddening 
distribution given by Neckel & Klare (1980). 



4 DUST PROPERTIES 

The wavelength dependence of polarization towards many 
Galactic directions follows the empirical relation (Serkowski 
et al. 1975; Coyne, Gehrels & Serkowski 1974; Wilking, 
Lebofsky, & Rieke 1982) 
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Figure 10. Upper panels: e versus radial distance of all stars from the center of Be 59 and e versus Pmax- Lower panels: ai versus radial 
distance of all stars from the center of Be 59 and ai versus Xmax- Cluster members are shown with open circles, foreground non- members 
with filled circles, stars in 3rd group with filled squares, late type stars with filled diamonds, other non-members with filled triangles, 
probable members with encircled filled circles and stars with uncertainty in their membership with encircled plus symbols. 



where Pa is the percentage polarization at wavelength A 
and Pmax is the peak polarization, occurring at wavelength 
Amax. The Xmax IS a functlon of the optical properties and 
characteristic particles size distribution of aligned grains 
(Serkowski, Mathewson & Ford, 1975; McMillan, 1978). The 
value of Pmax is determined by the column density, the 
chemical composition, size, shape, and degree and orienta- 
tion of the dust grains. The parameter K, an inverse mea- 
sure of the width of the polarization curve, was treated as a 
constant by Serkowski et al. (1975), who adopted a value of 
1.15 for all the stars. The Serkowski relation with 7^=1.15 
provides an adequate representation of the observations of 
interstellar polarization between wavelengths 0.36 and 1.0 
Hm. In one case Pmax and Xmax were obtained using the 
weighted nonlinear least square fit to the measured polar- 
ization by adopting; (1) K=1.15 for stars having data in 



V{RI)c passbands, or (2) A"=1.66 Xmax (Whittet et al. 
1992) for stars having data in B, V{RI)c passbands. Table 
[7] lists the Pmax, Xmax, o"i and e for 69 stars. The estimated 
values of the Pmax and Xmax using B, V{RI)c are listed 
in the second and third columns and those with V{RI)c 
passband data are in the seventh and eighth columns. We 
also computed the parameters ctJ^, (the unit weight error 
of the fit for each star) which quantifies the departure of 
the data from standard Serkowskis law and e, the disper- 
sion of the polarization angle for each star normalized by 
the average of the polarization angle errors (cf. Marraco, 
Vega, & Vrba 1993) . The estimated values of cti and e using 

t The values of cri for each star are computed using the expression 
(Tj = ^{r X / tp\)^ / (m — 2); where m is the number of colours and 

r\ = P\ Pmax eK\>[-K In^ {Xmax / X) ■ 
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Table 7. The Pmaxi ^7naxi f^l and e for the observed 69 stars 



37 stars with BV(/?/)c pass-band data 32 stars with Y{RI)c pass-band data 



Star ID* 




ax 


± 


e 


^TYlOiX 6 
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Star ID * 


Pmax 


ite 


^max i £ 




e 




(per cent) 












(per cent) 


(/i'm) 






(1) 




(2) 






(3) 


(4) 


(5) 


(6) 




(7) 






(8) 


(9) 


(10) 


3 


5.19 


-|- 


0. 


.08 


0.55± 0.02 


0, 


.85 


1.84 


47 


4.05 


-|- 


0.28 


0.69 


di 0.14 


0.57 


1.12 


15 


6.01 


-|- 


Q 


.12 


0.62± 0.03 


1, 


.14 


0.22 


59 


4.72 


-|- 


0.31 


0.68 


di 0.13 


0.33 


1.62 


16 


7.38 


-|- 





.09 


0.57± 0.02 


1, 


.67 


2.28 


62 


4.69 


-|- 


0.14 


0.58 


di 0.03 


3.29 


0.47 


16A 


5.57 


-|- 


0. 


.04 


0.52± 0.01 


3, 


.54 


3.08 


64 


6.02 


-|- 


0.20 


0.61 


di 0.04 


0.79 


0.45 


20 


4.36 


-|- 


0. 


.19 


0.55± 0.06 


0, 


.85 


0.33 


92 


6.00 


-|- 


1.10 


0.48 


lb 0.10 


0.83 


0.43 


23 


6.36 


-|- 


0. 


29 


0.53± 0.06 


0, 


.87 


1.63 
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6.25 


-|- 


0.24 


0.50 


di 0.02 


3.44 


0.67 


66 


7.18 


-|- 


0. 


.17 


0.56± 0.03 


0, 


.15 


0.79 


120 


3.91 


-|- 


0.16 


0.64 


-j- 0.07 


3.49 


2.13 


69 


6.67 


-|- 


0. 


32 


0.52± 0.06 


1, 


.12 


1.30 
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6.59 


-|- 


1.03 


0.45 


-j- 0.07 


1.68 


0.27 
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5.56 
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0. 


.12 
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0, 


.94 


3.99 
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0.28 


0.52 


-j- 0.02 
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0.61 
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0. 


.10 
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1, 


.82 


5.87 
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2.31 


-|- 


0.20 


0.52 


lb 0.06 


1.06 


2.37 
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7.08 


-|- 


0. 


.51 


0.49± 0.08 


0, 


.23 


1.45 
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2.08 


-|- 


0.14 


0.60 


lb 0.09 
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3.00 
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7.32 
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.17 


0.55± 0.03 
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.56 


2.86 
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5.54 
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0. 
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0, 


.52 
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5.41 


-|- 


Q 


.27 
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0, 


.69 


2.57 


240 
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-|- 


0.25 
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di 0.05 


1.38 


2.13 


239 


2.26 


-|- 


Q 


.22 
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1, 


.03 


0.46 
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0.57 
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0.97 
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.04 
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1, 
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4.98 
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0.33 


0.47 
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4.28 
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.29 
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5.40 
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0.20 


0.60 
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0.36 


0.59 
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6.02 
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.17 


0.49± 0.03 


2, 


.13 


1.00 


594 


1.95 


-|- 


0.10 


0.47 


di 0.03 


2.81 


0.75 


316 


2.99 


-|- 


0. 


.25 


0.46± 0.09 


1, 


.31 


1.10 
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8.73 


-|- 


0.39 


0.49 


di 0.02 


3.57 


0.41 


324 


7.10 


-|- 





.30 


0.51± 0.05 
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.35 


1.48 


742 


5.74 
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0.29 


0.61 
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0.94 


1.05 
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D.Ui 


It 
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Q AA 
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n 1 n 
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it U.Ui 
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1.48 
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.25 
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0.55 
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0, 


.39 


0.52 


1106 
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it 


0.47 


0.45 
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6.93 


± 


0, 


.13 


0.57± 0.03 


0, 


.74 


2.99 


1135 


8.58 


it 


0.55 


0.47 


it 0.03 


3.73 


1.86 


418 


7.82 


± 


0, 


.22 


0.55± 0.04 


0, 


.76 


3.00 


1190 


2.77 
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0.16 


0.55 
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1.05 
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.08 
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0, 


.75 


5.89 


1273 


5.74 
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0.57 
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0.35 
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± 


0, 


.06 
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2, 


.36 
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5.39 
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0.48 
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0.97 


1523 


2.40 


it 


0.15 


0.61 


it 0.09 
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0.49 


1067 


4.57 


± 


0, 


.20 


0.52± 0.06 


0, 


.99 
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* : From P08 



BV{RI)c passband data are listed in the fourth and fifth 
columns whereas those using V{RI)c passband data are in 
the ninth and tenth columns of Table [T] Further, the Pmax 
and Xmax for the 37 stars, which have data in the BV{RI)c 
passbands, have been calculated using the data of only three 
passbands V{RI)c and K—1.15. A comparison of the Pmax 
and Xmax values obtained using the three and four passband 
data is shown in Fig. [9l which manifests a good agreement. 

If the wavelength dependence of polarization is well rep- 
resented by the Serkowski law, ai should not be greater 
than 1.5 because of the weighting scheme. A higher value 
(> 1.5) could be indicative of intrinsic stellar polarization 
(Waldhausen et al. 1999, Feinstein et al. 2008). The polar- 
ization angle rotation with wavelength (e) also indicates the 
presence of an intrinsic polarization or a change of Xmax 
along the line of sight (Coyne 1974, Martin 1974). System- 
atic variations with wavelength in the position angle of the 
interstellar linear polarization of star light may also be in- 



dicative of multiple dust layers with different magnetic field 
orientations along the line of sight (Messinger et al. 1997). 
Following the above stated criteria, we consider stars with 
cri>1.5 and e>2.3 as probable candidates to have intrinsic 
polarization and/or polarization angle rotation. 

Fig. [To] plots radial distance of stars from the center of 
the cluster versus I (upper left panel), e versus Pmax (upper 
right panel) , the radial distance of stars versus ai (lower left 
panel) and cri versus Xmax (lower right panel). One can see 
that a significant number of stars show deviation from the 
normal distribution (14 stars have e> 2.3, 19 stars have ai 
> 1.5 and 5 stars have cti > 1.5 as well as e>2.3 (cf. Fig. 
1101 and Table [7]). It is interesting to mention that the ma- 
jority of foreground stars (12 out of 14) have e<2.3 and all 
the stars having e>2.3 are located within the cluster region, 
whereas the majority of the stars (10 stars) having ai > 
1.5 are located within the cluster region. About 40 per cent 
(28 out of 69) stars of the sample show the signatures of 
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Figure 11. (J — H') versus (H — K) colour-colour diagram for 
all the observed stars towards Be 59. The data are taken from 
the Two-Micron All-Sky Survey (2MASS) Point Source Catalog 
(Cutri et al. 2003). 2MASS data have been converted into the 
California Institute of Technology system using the relations pro- 
vided by Carpenter (2001). The theoretical tracks for dwarfs and 
giants are drawn (Bessell & Brett 1988). Reddening vectors are 
also drawn (Cohen et al. 1981). The symbols are same as in that 
of the Fig. 11 



either intrinsic polarization or rotation in their polarization 
angles. Ten stars (#3, 16A, 133, 169, 195, 234, 367, 454, 
689 and 1106) are found to be located in the near infrared 
excess zone (see Fig. [TT|) . Six of them (# 16A, 133, 195, 
454, 689 and 1106) show intrinsic polarization, with tJi > 
1.5. The upper right panel of Fig. 1 101 indicates that a signif- 
icant number of stars (12 stars) with Pmax > 4.5 per cent 
have e > 2.3. The star #247 (BD-|-66° 1673) is located to- 
wards the north western edge of the cluster. This star was 
classified previously with objective prism spectra as 09— BO 
(MacConnell 1968; Walker 1965), but is reclassified later as 

05 V((f))n (Majaess et al. 2008), making it the hottest star 
in Cep OB4. The star's high temperature drives mass loss 
via strong stellar winds (Yang and Fukui 1992; Gahm et al. 
2006). Our polarimetric results indicate an intrinsic nature 
of polarization as it has cri = 1.83 and e=4.98. 

Another criterion to detect intrinsic stellar polarization 
is based on \max- A star having \max much lower than the 
average value of the ISM (0.545 pim; Serkowski, Mathewson 

6 Ford 1975) is considered as a candidate to have an intrinsic 
component of polarization (Orsatti, Vega, & Marraco 1998) . 
In the present study only one star #215 has been found to 
have a much lower value of Amaa;=0.36±0.05^m. This star 
has ai=5.5. 



4.1 Extinction law 

To study the nature of the extinction law in the cluster re- 
gion we used the TCDs as described by Pandey et al. (2000, 
2003), in the form of (V^ — A) vs. [B — V), where A is one of 
the wavelengths of the broadband filters R, I, J, H, K or L, 
to separate the influence of the normal extinction produced 
by the diffuse ISM from that of the abnormal extinction aris- 
ing within regions having a peculiar distribution of dust sizes 
(cf. Chini & Wargau 1990, Pandey et al. 2000). The (V - A) 
vs. {B — V) TCDs for the cluster region are shown in Fig. 
1121 The distribution of field stars follows a normal reddening 
law for the foreground inter-stellar matter. The slopes of the 
distributions for cluster members, niciuster, are found to be 
1.40±0.08, 2.62±0.15, 3.08±0.20, and 3.31±0.23 for (V-I), 
{V - J), {V - H), {V - K) versus (B - V) TCDs respec- 
tively. The ratios ^^j^Zy) ^^'^ ratio of total-to-selective 
extinction in the cluster region, Rduster, is derived using 
the procedure given by Pandey et al. (2003). Assuming the 
value of Rv for the diffuse foreground ISM as 3.1, the ratios 
E(B-v) yisld Rduster = 4.0±0.1, which indicates an anoma- 
lous reddening law. P08 have also estimated an anomalous 
reddening law in the cluster region with Rduster = 3.7±0.3, 
with a normal reddening law for the foreground diffuse ISM. 
In the central region of Be 59, MacConnell (1968) also found 
evidence for a large value (3.4—3.7) of Rv- Several studies 
have already pointed out an anomalous reddening law with 
a high Rv value in the vicinity of star-forming regions (see 
e.g. Pandey et al. 2003 and references therein), however for 
the Galactic diffuse ISM a normal value of Rv = 3.1 is 
well accepted. The higher than the normal values of Rv has 
been attributed to the presence of larger dust grains. There 
is evidence that within dark clouds accretion of ice mantles 
on grains and coagulation due to colliding grains change 
the size distribution towards larger particles. On the other 
hand, in star-forming regions, radiation from massive stars 
may evaporate ice mantles resulting in small particles. Here 
it is interesting to mention that Okada et al. (2003) , on the 
basis of the [Sill] 35 iim to [Nil] 122 /xm ratio, suggested that 
efficient dust destruction is occurring in the ionized region of 
Be 59. Chini & Kruegel (1983) and Chini & Wargau (1990) 
have shown that either larger or smaller grains may increase 
the ratio of total-to-selective extinction. 

The weighted mean value of \max for the cluster re- 
gion was estimated to be 0.538±0.004^m. The mean value, 
within error, is comparable to the value measured in the 
general ISM (0.545 ^m, Serkowski et al. 1975). Using the 
relation Rv = (5.6±0.3)A„a^ (Whittet & Van Breda 1978), 
the value of Rv, the total-to-selective extinction, comes out 
to be 3.01±0.16, which is in agreement with the average 
value [Rv ~ 3.1) for the Milky Way, but is in contradic- 
tion with the result obtained from {V - X)/{B - V) TCDs. 
The mean value of Xmax for foreground stars is estimated as 
0.498±0.017Atm, which yields the Rv value as 2.79±0.18 for 
the foreground diffuse ISM. This indicates a smaller value 
of Rv for the diffuse ISM towards the direction of the Be 
59 {I — 118.2°). There is much evidence in the literature 
that indicates significant variations in the properties of inter- 
stellar extinction along various Galactic directions. Whittet 
(1977) reported that the value of Rv in the Galactic plane 
can be represented by a sinusoidal function of the form Rv ~ 
3.08-1-0.17 sin(Z + 175°), which indicates a minimum value 
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Figure 12. {V — I), {V — J), {V — H), {V — K) versus (B — V) two-colour diagrams for the observed stars towards Be 59. The first 
group of foreground stars are shown with filled circles, which follow the normal reddening by the general diffuse ISM (see Table 3 of 
Pandey et al. 2003). The cluster members (open circles) show anomalous reddening as their slopes are significantly different from those 
of foreground stars. The fitted slopes for cluster members at each wavelength combination are over laid on the figure. 



of Rv at I 95°. The above relation suggests a value of 
Rv ~ 2.9 towards the direction of Be 59 (i.e. I = 118.2°). 
The study by Geminale & Popowski (2004) also suggests a 
lower value (^ 2.9) towards Galactic longitude I ~ 120°. 
Thus the present estimation of Xmax, and consequently the 
value of Rv for the diffuse ISM is in agreement with the 
values reported in the literature. Hence we conclude that 
the Rv for the intra-cluster matter could be higher in com- 
parison to that for the general diffuse matter towards the 
direction of Be 59. Here it is worthwhile to note that there 
is much evidence of variation of grain size-distribution to- 
wards the direction of Be 59. 



4.2 POLARIZATION EFFICIENCY 

The ratio of Pmax / E{B — V) is known to be a measure 
of the polarization efficiency of the ISM and it depends 
mainly on the grain alignment efficiency, the magnetic field 
strength and the amount of depolarization due to the radi- 
ation traversing clouds with different magnetic field direc- 
tions. When the light passes through multiple dust layers. 



the resultant polarization may increase or decrease depend- 
ing upon the orientation of the magnetic field in each dust 
layer (see for example, Feinstein et al. 2003b, Orsatti et al. 
2003, Martinez et al. 2004, Vergne et al. 2007, 2010 and Ell) 
The observed polarization and extinction data towards a 
particular direction of the Galaxy provides important inputs 
to test the models dealing with the extinction and alignment 
of the grains. Figure displays the polarization efficiency 
diagram for the observed stars. The symbols are the same as 
in Fig. 1101 It is well known that for the diffuse ISM the po- 
larization efficiency can not exceed the empirical upper limit 
given by, P™ax=3 Av ~ 3 Rv E{B-V) ~ 9.3 E{B-V) per 
cent (assuming i?v=3.1, Hiltner & Johnson 1956; Serkowski 
et al. 1975) and the same is shown by a continuous line in 
Fig. [I3l For the average ISM, Serkowski et al. (1975) have 
found that the polarization efficiency of the ISM follows the 
mean relation Pmax — 5 E{B — V), which is shown by a 
dashed line. The recent estimate of the average polariza- 
tion efficiency for the general diffuse ISM by Fosalba et al. 
(2002), which is valid for E{B - V) < 1.0 mag, is shown 
with a dash-dotted line. 

Figure [13] shows that the foreground stars (filled circles) 
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Figure 13. Pmax versus E{B — V) for the stars with available reddening (cf. Table (Tjl. The symbols are same as that of the Fig. 1101 
The stars with two possible reddening values are connected by thin dash-dotted line. The solid line represents the empirical upper limit 
relation for the polarization efHciency (by assuming _Rv=3.1) of Pmax = 9.3 X E{B — V) (Serkowski et al. 1975). The dashed line 
represents the relation Pmax = 5 X E(B — V) (Serkowski et al. 1975) and the dash-dotted line represents the relation Pmax — 3.5 X 
E(B - y)0-* by Fosalba et al. (2002). 



are distributed along the dashed hne, which suggests that 
the dust grains in the dust layer located at < 470 pc have 
polarization efHciency comparable to the average polariza- 
tion efficiency (~5 per cent per mag) of the diffuse ISM. The 
stars located at ~ 700 pc are shown with filled squares. The 
colours of these stars indicate two values of E{B — V) (see 
Table We have used both values of E{B - V) for these 
stars and those data are connected with thin dashed-dotted 
Unes. In general it seems that the dust layer located at 500 
< d < 800 pc also has an average polarization efficiency (i.e., 
5 per cent per mag). The majority of the cluster members 
(open circles) are distributed below the dashed line thereby 
indicating that the intra-cluster medium exhibits less polar- 
ization efficiency than the mean value for the diffuse ISM (~ 
5 per cent per mag). The large dispersion in Pmax (4—8 per 



cent) for cluster members is compatible with the differential 
reddening within the cluster (~ 1.4—1.7 mag). 

The net polarization due to the intra-cluster medium 
is estimated to be ~2.2 per cent (cf. Table E)). The differ- 
ential E{B — V) due to the intra-cluster medium is ~0.3 
mag. Thus the net polarization efficiency due to the intra- 
cluster medium comes out to be ~ 7.3 per cent per mag, 
which is higher than that due to the diffuse ISM. The small 
dispersion in the mean value of 6v (4°, Table [5]) of cluster 
members also indicates a better alignment of dust grains in 
the intra-cluster medium. As discussed in Sec 13.31 for the 
foreground stars of the first group and third group as well 
as the cluster members, the mean polarization angle changes 
with increasing distance systematically which may lead to 
the depolarization effect in the case of radiation from clus- 
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ter members. Hence the less polarization efficiency of the 
intra-cluster medium, as seen in Fig. 1131 could be because 
of different alignment of dust grains in the foreground dust 
layers. Similar kind of depolarization effect has been ob- 
served towards Trumpler 27 by Feinstein et al. (2000), Hogg 
22 and NGC 6204 by Martinez, Vergne and Feinstein (2004) 
and towards NGC 6124 by Vergne et al. (2010). Here it is in- 
teresting to mention that stars # 418, 441 and 461, located 
at the northern periphery of the cluster have the highest po- 
larization efficiency among the cluster members. These stars 
have polarization efficiencies greater than the average value 
(5 per cent per mag) for the diffuse ISM. Four stars # 689, 
828, 1106 and 1135 are located outside the boundary of the 
cluster and towards north-east of the cluster center. These 
stars have the highest polarization {Pv ~ 7.6—8.7 per cent) 
and polarization efficiency greater than 5 per cent per mag. 
The E{B — V) values of these stars range ~1. 20— 1.50 mag. 
Interestingly, all these four stars have their tJi > 1.5 (cf. Fig. 
llOf) . thereby indicating presence of intrinsic polarization. 



5 SPATIAL VARIATION OF E{B - V), Pv, Ov 

AND \max 

Fig. [T4ja) shows the spatial distribution of E{B — V) as 
a function of radial distance from the center of the clus- 
ter. The distribution reveals a lack of reddening material 
near the center of the cluster. The E{B — V) increases away 
from the center. The E{B — V) values reaches its maximum 
at ~ 4' (~ 1.2 pc) and remains constant up to ^ 6' (1.8 
pc). The E{B — V) value decreases for radial distances > 
6'. The spatial distribution suggests that the density of the 
reddening material is high at ~ 1.2—1.6 pc in comparison to 
that near the cluster center. On the basis of one-dimensional 
raster-scan observations of Shl71 (Be 59) region, Okada et al 
(2003) have studied the the spatial distribution of line inten- 
sities, representing the lowly-ionized gas, the highly ionized 
gas and photo dissociation region (gas phases) as a function 
of radial distance from the cluster. The lines representing 
the neutral region (e.g., [Ol] 63/im, 146/im; [Cll] 158/im; 
[Sill] 35/xm), lowly ionized region ([Nil] 122/xm) and highly 
ionized region ([Olll] 52/im, 88/xm; [Nlll] 57/im) show max- 
ima at ~ 1 pc and 4 pc indicating the presence of high 
density gas at these positions. Their figure 6, showing 0.61 
GHz continuum flux on the line along the raster-scan obser- 
vations, also indicated that the flux increases systematically 
with radial distance from the cluster center and reaches the 
maximum at ^ 1.2 pc and remains constant up to ^ 1.7 
pc. The flux shows a decreasing trend for radial distance > 
1.7 pc. The spatial distribution of 0.61 GHz continuum flux 
agrees nicely with the E{B — V) distribution as shown in 
Fig. ll4f aV Thus, it is clear that the density of the gas in- 
creases with the radial distance from the cluster center with 
the maximum at ~ 1 pc. 

The spatial variation of Pv as a function of radial dis- 
tance from the cluster center (Fig. ll4r b)) shows a systematic 
increasing trend with the radial distance up to ~ 5'; how- 
ever for radial distance >5', the distribution of Pv shows no 
trend but a significant scatter. The variation of Pv agrees 
well with the distribution of gas as discussed above, indi- 
cating that the average polarization efficiency in the core as 
well as in the corona of the cluster is the same. Fig. I14f c) 



shows the distribution of dv as a function of radial distance. 
The 6v distribution also shows a systematic change with the 
radial distance, in the same sense that the average 9v value 
(~ 110°) at the center decreases to ~100° at ~5'. For ra- 
dial distances >5', the distribution of Ov shows a scattered 
distribution around 102°. On the basis of the 9v distribu- 
tion we can conclude that the magnetic field orientation in 
the core of the cluster is significantly different from that 
in the coronal region, which is comparable to the magnetic 
field orientation of the third group of the stars (~ 96°). We 
presume this could be either due to the molecular cloud at 
the center being perturbed during the cloud collapse or to 
strong stellar winds/supernova explosion. 

Fig. I14f dl shows the spatial variation of Xmax obtained 
for identified members of the cluster, which manifests that 
the Xmax is higher near the center of the cluster. The av- 
erage value of Xmax decreases up to ~ 5' with increasing 
radial distance from the cluster center. The distribution of 
Xmax suggests that the dust grain size near the center of 
the cluster is relatively higher in comparison to that in the 
coronal region. Okada et al (2003), on the basis of the [Sill] 
35/im to [Nil] 122/xm ratio, suggested that the efficient dust 
destruction is occurring in the ionized region. It is possible 
that the smaller dust grains might have been evacuated from 
the central region of the cluster due to strong stellar winds 
leaving relatively large dust grains in the central region of 
the cluster. 



6 CONCLUSIONS 

In the present study we have carried out polarimetric 
observations towards the direction of young open cluster 
Be 59 in the B,V,{R,I)c bands. The aim of the study 
was to investigate the properties of dust grains in the ISM 
towards the direction of Be 59 as well as the properties of 
intra-cluster dust. 

Following are the main conclusions of the present study: 

The distribution of Pv and dv suggests three dust layers 
towards the direction of Be 59 at ~300 pc, ~500 pc and 
~700 pc respectively. The total polarization due to these 
dust layers is found to be ~0.2— 1.0, ~1. 0—3.0 and ~ 5.5 
per cent respectively. The magnetic field orientation in these 
dust layers is different from each other. The magnetic field 
orientation of the first dust layer ('^ 82°) is rather similar 
to that of the GP (86°). 

We have further shown that the polarization measure- 
ments in combination with the {U — B) — {B — V) colour- 
colour diagram provide a better estimation of the cluster 
members. The polarization measurements of the identified 
cluster members reveal that the net polarization due to 
the intra-cluster medium is estimated to be ~2.2 per cent. 
About 40 per cent cluster members show the signatures of 
either intrinsic polarization or rotation in their polarization 
angles. 

The TCDs of the identified cluster members reveal 
an anomalous reddening law in the cluster region. The 
weighted mean values of Xmax for the cluster region and 
for the foreground stars are estimated as 0.54±0.01/xm and 
0.50±0.02/im. The above estimate of Xmax for the fore- 
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Figure 14. Radial variation of E{B — V), Py, 9v and Xmax for only the cluster members. In the top panel, only the cluster members 
with E{B — V) are used. The estimated error (0.05 mag) in E{B — V) is also plotted for a reference. 



Polarimetric study of Berkeley 59 21 



ground ISM indicates relatively smaller dust grain sizes, 
consequently smaller value of total-to-selective absorption 
_Rv=2.79±0.18 in comparison to the normal value for the 
diffuse ISM {Rv~S.l). Thus mean Xmax value for the clus- 
ter region (0.54±0.01/im) suggests a relatively larger grain 
size in the cluster region in comparison to those in the gen- 
eral diffuse ISM towards the cluster region. 

The foreground dust layers have polarization efHciency 
comparable to the average polarization efficiency (~ 5 per 
cent per mag) of the diffuse ISM, whereas the majority of the 
cluster members indicate a smaller polarization efficiency for 
the intra-cluster medium. It indicates that the star light of 
the cluster members might have been depolarized because of 
the non-uniform alignment of dust grains in the foreground 
dust layers and in the intra-cluster medium. 

The spatial distribution of E{B — V) in the cluster re- 
gion shows an increasing trend with the radial distance. 
The polarization is also found to be systematically increas- 
ing with the radial distance from the cluster center. Both 
E{B — V) and Py values reach a maximum value at ~4'- 
5' {'-^ 1.2—1.5 pc). The 9v as well as Xmax for the cluster 
members are found to decrease systematically with increas- 
ing radial distance. 
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